Unusual Brownian motion of photons in open absorbing media 
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Very recent experiments have discovered that localized light in strongly absorbing media displays 
intriguing diffusive phenomena. Here we develop a first-principles theory of light propagation in 
open media with arbitrary absorption strength and sample length. We show analytically that 
photons in localized open absorbing media exhibit unusual Brownian motion. Specifically, wave 
transport follows the diffusion equation with the diffusion coefficient exhibiting spatial resolution. 
Most strikingly, despite that the system is controlled by two parameters - the ratio of the localization 
(absorption) length to the sample length - the spatially resolved diffusion coefficient displays novel 
single parameter scaling: it depends on the space via the returning probability. Our analytic 
predictions for this diffusion coefficient are confirmed by numerical simulations. In the strong 
absorption limit they agree well with the experimental results. 



PACS numbers: 42.25.Dd,71.23.An 

Many complex systems exhibit fundamentally differ- 
ent physics in various scales. A canonical example, thor- 
oughly studied by Einstein in one of his revolutionary 
papers in 1905 [T, is the motion of pollens suspended 
in water (the Brownian motion). In short time scales, 
pollens undergo bombardment by water molecules and 
exhibits random motion. The resulting complexity of mi- 
croscopic dynamics notwithstanding, in large time scales 
emerges substantially simpler physics: the evolution of 
the pollen concentration profile fully characterizes the 
dynamics and is governed by a diffusion equation. Such 
emergent diffusive phenomenon is universal, and micro- 
scopic dynamics only affect the diffusivity namely the 
value of the diffusion coefficient. 

The classical electromagnetic wave in infinite random 
media is the very system where universal macrosocpic dif- 
fusive phenomena emerge from complicated microscopic 
dynamics. In short (time) scales, the system's physics is 
dominated by light scattering off random dielectric field 
described by the Maxwell equation. In scales larger than 
the elastic mean free time, as noticed by astrophysicists 
long time ago [5], light propagation may be viewed as 
random motion of Brownian particles, the photons. Yet, 
light displays the particle- wave duality, and the wave na- 
ture has far-reaching consequences. Specifically, wave in- 
terference suppresses the diffusivity [3]. In the presence 
of strong dielectric fiuctuations, the (bare) diffusion co- 
efficient is strongly renormalized and vanishes in a much 
larger time scale. This is a hallmark of Anderson lo- 
calization of light (see Ref. [4; for a review), and the 
corresponding time scale is called the localization time. 

In realistic experiments or photonic devices the 
medium is open: photons escape the medium from the 
interfaces. Therefore, diffusivity of waves cannot be fully 
suppressed even disorders are strong. A fundamental 
question thereby is: do localized waves in open media 



exhibit certain macroscopic diffusive phenomena? This 
is a long-standing issue (see Ref. |5j for a recent review) . 
Conceptually, a prominent difficulty arises from the fol- 
lowing. On the one hand long-time transport of localized 
waves in open media is dominated by rare long-lived, high 
transmission (resonant) states [6l[7], while on the other 
hand macroscopic diffusion of waves - if exists - char- 
acterizes statistical behavior of wave propagation: these 
are two seemingly contradictory facets. In fact, experi- 
ments and numerical simulations [SVllOj have shown that 
due to the failure of capturing resonant states the pre- 
vailing macroscopic diffusion model developed in Ref. [11] 
cannot describe transport of strongly localized waves in 
open media. This had inspired the opinion that the reso- 
nant state and macroscopic diffusion of photons are two 
incompatible concepts. 

Contrary to this common belief, the results achieved 
recently |H] show that the rare resonant states do not 
wash out macroscopic diffusion of photons. Rather, they 
are reconciled in the way that the spatial resolution of the 
diffusion coefficient exhibits a novel scaling. Specifically, 
the diffusion coefficient depends on the space via the re- 
turning probability. (One may show - both analytically 
and numerically [9, - that this novel scaling is completely 
missed in the prevailing model [llj due to the unjustified 
self-consistent treatments employed there.) This novel 
scaling was analytically found by using a first-principles 
theory [51 112) and has been fully confirmed by numerical 
simulations [S]. 

So far we have focused on passive media. Yet, absorp- 
tion inevitably exists in realistic environments |13| . In 
open media the strong interplay between absorption and 
localization has been experimentally shown to cause sur- 
prising phenomena |14j . and may even find practical ap- 
plications in optical devices [15] . Very recently, both real- 
istic and numerical experiments have observed that local- 



ized waves exhibit unusual macroscopic diffusion in open 
absorbing media |16l 117) . However, these experiments 
(both realistic and numerical) focus on short samples and 
strong absorption. For such systems resonant states are 
unimportant. In the opposite limit, where samples are 
sufficiently long and absorption is weak, resonant states 
dominate transport of localized waves. As these states 
are very sensitive to absorption due to long lifetimes [18] , 
a very fundamental issue arises: how do absorption in- 
terplay with macroscopic diffusion of localized waves in 
open media? The purpose of this Letter is to present 
the first systematic first-principles study of this issue for 
arbitrary sample length and absorption strength. 

Main results and experimental relevance. — Specif- 
ically, we consider propagation of light in quasi-one- 
dimensional open (uniformly) absorbing media of length 
L. This system is controlled by two dimensionless pa- 
rameters, £,/L and ^a/L, where ^ (^q) is the localization 
(absorption) length. We find that - even in localized 
samples {£,/L <C 1) - photons exhibit macroscopic diffu- 
sion irrespective of the absorption rate 7. More precisely, 
in the presence of a steady monochromatic light source 
located at x' the wave intensity (namely energy density) 
profile, y{x,x'), obeys 



?a/e, 



{7 - d^D{x)d,} y{x, x') = 6{x - x'). 



(1) 



It resembles the normal diffusion equation, but, the diffu- 
sion coefficient exhibits a number of anomalies. First of 
all, it has a spatial resolution, D{x), decreasing monoton- 
ically from the air-medium interface to the sample mid- 
point and varying over several orders. Most strikingly, 
despite that the system is controlled by two parameters, 
D{x) exhibits a novel single parameter scaling. 



D{x)/DiO)=DooiXix)) 



(2) 



similar to the passive medium case. Here, the scaling fac- 
tor A(2;) is proportional to the (static) returning proba- 
bility which depends on x/L as well as ^/L and £,a/L, and 
the scaling function Doo(A) is the same as that of pas- 
sive media. (The explicit form of both A(a;) and Doo{X) 
will be given below.) Bearing these anomalies and the 
analogy to the normal diffusion equation, Eq. (nj) de- 
scribes 'unusual Brownian motion' of photons. Further- 
more, from Eq. (pi) we are able to predict analytically 
the spatial resolution of the diffusion coefficient which 
is confirmed by numerical experiments. We stress that 
our analytic results are very general: they are valid for 
arbitrary sample length and absorption strength and, ad- 
ditionally, for both time-reversal (orthogonal symmetry) 
and broken time-reversal (unitary symmetry) systems. 

In the case of strong absorption, our analytic results 
are simplified and fully agree with the recent experimen- 
tal findings shown in Refs. [16l[T7j. Specifically, if absorp- 
tion is strong such that ^a/L ^ 1, we find that D{x) ex- 
hibits a plateau in the regime min(a;, L — x) ^ ^a- More- 
over, the plateau value is determined by single parameter. 



DiL/2)/D{0)^D^{U2O 



(3) 



(cf. the first line of Table [l] and [h]) . We stress that the 
plateau as well as Eq. (Is]) exists only ii £_a/L <S^ 1. 

Origin of unusual Brownian motion. — We begin 
with a physical explanation of the main results. Since 
closer to the interface more easily escape photons from 
the medium (cf. dashed lines of Fig. II]), the returning 
probability, yf){x,x), is inhomogeneous in space (i.e., de- 
pending on X as well as the parameters ^/L and ^a/L). 
As such, wave interference effects are inhomogeneous 
also. Near the interface (deep inside the sample) they 
are weak (strong) . For diffusive samples {£,/L :^ 1 where 
£_ ^ ttvDq [19] with v the photon density of states), 
yfi{x,x) gives rise to an inhomogeneous (one-loop) weak 
localization correction ^ 0{yo{x,x)/{TTv)) ^ 0{\{x)). 




FIG. 1: Examples of wave interference between two opti- 
cal paths (red and blue) leading to the novel scaling behav- 
ior. The paths trace the three loops with different orders 
(B-5>A->C for the red and C->B->A for the blue). Notice 
that these loops may not be complete (e.g., dashed parts of 
A and C) due to photon leakage. 

For localized samples (^/L ^ 1), as waves penetrate 
deeply into the sample photons easily return to its de- 
parture point. This causes more complicated wave in- 
terference essential to strong localization physics. For 
example, two optical paths may take the same n{— 3 in 
Fig. IlJ loops as their routes and trace them with differ- 
ent orders. Therefore, they constructively interfere with 
each other. This leads to a weak localization correction 
~ C'(A"(a;)). Therefore, the spatial resolution D{x) de- 
pends on x/L (as well as ^/L and £,a/L) via the factor 
A(x), justifying Eq. ^. 

Importantly, the optical paths in Fig.[T]may propagate 
in the same direction - clockwise or counterclockwise - 
during tracing each loop, suggesting that interference pic- 
ture beyond what leads to usual one-loop weak localiza- 
tion does not necessarily require time-reversal symmetry. 
This is crucial to localization physics of systems with uni- 
tary symmetry (which is completely beyond the reach of 
the self-consistent theory [Hj.) Therefore, the novel sin- 
gle parameter scaling is an intrinsic wave phenomenon 
unrelated to time-reversal symmetry. 

First-principles theory. — The derivations are largely 
parallel to those of passive media [HI [H]. Therefore, 
we shall outline the key steps below with an emphasis 
on the key differences, and refer the readers to Ref. |5] 



for technical details. We first introduce the microscopic 
formalism valid for arbitrary dimensions. Consider the 
point-like source located at r' inside the medium with 
the spectral decomposition Ji^(r') (w being the wave an- 
gular frequency). The time-integrated wave intensity 
- upon disorder averaging - is given by [2U] 7(r) = 
/ |^J^(r,r')|J„(r')p. Here the (static) spatial correla- 
tion function 3^(r, r') = (G^^ (r, r') G^^ (r', r)), with (• • • ) 
standing for the disorder average and the advanced (re- 
tarded) Green function G^ {G^) defined as 

[V2 + c.2(l + e{v) ± ze")]Gfi^(r, r') = ^(r - r'), (4) 

where the dielectric fiuctuation, e(r), is gaussian, e" > 
causes uniform absorption, and the light velocity in the 
air is set to unity. 

Then, it is a standard procedure to cast the above spa- 
tial correlation function in terms of the functional inte- 
gral over a supermatrix field (5(r), 

y{r,r') = ["^y J D[Q]str{A+Q{r)A^Q{r'))e-''^'^l{5) 

Here A± are some constant supermatrices and 'str' the 
supertrace [3 [12] . The action. 



F[Q] = ^ I dvsti{D^{VQf - 27AQ}, 



(6) 



differs from the passive one in the second term ac- 
counting for the absorption. Notice that the (bare) 
diffusion coefficient Dq = D{0) and for e" ^ 1 
the absorption rate 7 = cje" 21J. The supermatrix 
Q = T~^KT where A is the so-called metallic saddle 
point and T takes the value from the coset space of 
UOSP(2,2\A)/UOSP{2\2) ® UOSP[2\2) for orthogonal 
symmetry and of C/(l, 1|2)/C/(1|1) ® U{l\l) for unitary 
symmetry. It is very important that if the air-medium 
interface is transparent, the supermatrix field is fixed to 
be A at the interface. This reflects that on the interface 
there is no photon accumulation [5J [HI [T^ . 

In quasi one dimension the microscopic formalism is 
simplified. Specifically, the supermatrix Q field is ho- 
mogeneous in the transverse plane. F[Q] reduces to 
^ Jq dxstY {DQ{dxQY ~ 27A(3}, [v includes a factor of 
the cross section area.) and the boundary condition to 



Q{x = 0) = Q{x = L) = A. 



(7) 



As a result, y depend only on the longitudinal coordi- 
nate, X. Then, we follow the procedures of Refs. [HIH] to 
calculate y{x,x') exphcitly. We find that y{x,x') satis- 
fies Eq. ^ with the boundary condition: y{x — 0,x') — 
y{x = L, x') — 0. Furthermore, we find that the diffusion 
coefficient has a spatial resolution D{x), and the latter is 
a functional of the factor X{x) = yo{x,x' — x)I{'kv). 
This justifies Eq. ([2|. The correlator yo(x,x') solves 
the normal diffusion equation: (7 — Doi9^)3^o(3;,a;') = 



5{x~x') implemented by the boundary condition yo{x = 
0, x') — y^lx = L, x') = 0. Solving this equation we ob- 
tain 

X / X _ Ca COsh(L/^a) - COsh((L - 2x)/ia) ,^. 



2^ 



ML/ia) 



where the diffusive absorption length S,a = ^-Dq/t- 

For A ^ 1 we find that the perturbative expansion of 
the scaling function £'oo(A) is identical to that at 7 = 0. 
This implies that Z)oo(A) is the same as the one at 7 = 0. 
On the other hand, at 7 = the scaling function has 
been found analytically and fully confirmed by numerical 
experiments [S] , which has the following asymptotic form: 



DooW 



1 + CiA + C2A2 + --- , 

e-^, A>1. 



A<1, 



(9) 



Notice that the coefficients Ci depend on the system's 
symmetry. In particular, for orthogonal symmetry ci < 
while for unitary symmetry ci =0, C2 < 0. 

Eqs. Q, ([2|, ([8]) and ^ constitute a complete de- 
scription of unusual Brownian motion of photons in open 
absorbing media. In essence, it is a macroscopic phe- 
nomenon emerging from the microscopic Helmholtz equa- 
tion H in large time scales. The present macroscopic 
theory differs from the previous one for passive media 
O IH] in the absorption term and that the spatial res- 
olution D{x) depends also on a new parameter £,a/L. 
Below we provide numerical evidence of this intriguing 
phenomenon. 




FIG. 2: In localized open absorbing media the spatial resolu- 
tion of the diffusion coefficient exhibits novel single parameter 
scaling. Analytic predictions for D{x) (solid lines) are in good 
agreement with the simulation results (circles). From bottom 
to top the absorption strength increases (see the text for the 
value of e" as well as other parameters). 

Evidence from numerical experiments. — We perform 
numerical simulations of the (microscopic) Helmholtz 
equation. To this end we prepare a randomly layered 



medium where the analytic results remain valid. The rel- 
ative permittivity in each layer is a random number with 
a real part uniformly distributed in the interval [0.3, 1.7]. 
The imaginary part, e", is assumed to be a constant. The 
random medium is embedded in the air background so 
that there exist no internal reflections when e" vanishes. 
We launch a plane wave of frequency lo — 1.65a^^ into 
the system, where the sample length L = 200a with a 
the layer thickness, and use the standard transfer matrix 
method to calculate the wave intensity profile, Iii{x), for 
each dielectric disorder configuration /3. For a given value 
of e", we calculate the ensemble-averaged wave intensity 
profile, I{x) = {Ij3{x)), of 5 000 000 disorder configura- 
tions. Then, by presuming Eq. (fl]), which gives 



TABLE I: Behavior of D{x) at ^ < ^„ < L 



Dix) 



-J l^ I{x')dx' + D{0)d,\.,^Ll{x) 
dxl(x) 



(10) 



we obtain numerical results of D{x) from the measure- 
ments oi I{x). 

The results of simulations are shown in Fig. [2] for e" = 
0, 1 X lO"'', 5 X 10-"*, 1 X 10-3, and 2 x 10"^ (from bot- 
tom to top). First of all, in the passive limit (e" = 0, the 
bottom curve) we recover the result of Ref. [1] for deeply 
localized (Indeed, in this limit L/^ = 200a/22a « 9 is 
large.) samples. The effects of absorption become signif- 
icant when ia/L - 0(1), i.e., e" « ^/{L'^lo) w 3 x 10"^. 
Having this estimation in mind, wc may adjust the value 
of e" to systematically explore the effects of absorption. 
The results are compared with the analytic predictions 
obtained from the novel single parameter scaling theory 
namely Eqs. ([2|, Q and (|9|. As shown in Fig. [2] they 
are in good agreement. Now we turn to analytic dis- 
cussions of the behavior of D{x) based on the developed 
first-principles theory. 

Behavior of D(x) in semi-infinite samples. — In this 
case L ^^ (X), the scaling factor (Is]) is simplified to 
\{x) = ^{^ ~ e"^^/^"). If the absorption is weak such 
that ^/^a < 1, then A < 1 for x < ^ and A > 1 for a; > ^. 
Applying the first and second line of Eq. ^ , respectively, 
we find distinct behavior of D{x)/Do in various regimes, 
which is summarized in Table HI The first line indicates 
that D{x) reaches an exponentially small residual value 
at a; > Ca/2. This reflects that deep inside the medium 
the ordinary Brownian motion of photons is restored, al- 
beit with dramatically smaller diffusion coefficient. Qual- 
itatively, the plateau arises from that due to absorption 
photons could not propagate for a distance larger than 
^a and as such, photons deep inside the medium could 
not 'see' the air-medium interface. Indeed, similar re- 
sults were achieved by using the scaling theory of light 
localization in infinite absorbing media long time ago [3] . 
The second and third lines indicate an intermediate scale, 
y/£,a£., between ^ and £,. For x smaller than this scale the 
local diffusion coefficient behaves essentially the same as 
that of (semi-infinite) passive media (see the third and 



regime 


Dix)/D{0) 


X > ia/2 


g-€a/{2e) 


V^a^ < 2; < ^a/2 


f,^^/iiaOg-^/i 


^<x<^^M 


e-^/« 


a;<C 


1+Ci(x/C)+C2(x/C)"+--- 



fourth lines), while the deviation starts at x ~ V^aS.- A 
physical explanation for this deviation may be as follows. 
Consider light incident from the interface. It effectively 
penetrates into the medium with a depth f^. As such, 
the Lyapunov exponent (inverse localization length), 7, 



fiuctuates, following a distribution ~ e 



U-/-i-'f 



Averaging e '^'^ with respect to this distribution we re- 
cover the result given in the second line. 

If absorption is strong such that S,/£,a ^ 1, then A <C 1 
irrespective of x. By using the first line of Eq. ([9]) we find 
the results summarized in Table ITll 

Behavior of D(x) in finite localized samples. — These 
samples {£,/L <C 1) have an essential difference from semi- 
infinite samples in the existence of resonant states most 
of which reside near the sample center O [7] . As shown 
in Ref. [5], they play decisive roles in establishing the 
highly unconventional macroscopic diffusion of localized 
waves. Due to long lifetime the resonant state is very sen- 
sitive to absorption. Below wc shall discuss separately 
the cases of weak {S,a/L » 1) and strong {£,a/L <^ 1) 
absorption. In the former case, the scaling factor (Is]) is 
simplified to X{x) « ^^^^[1 - i^^^)]. For A(a;) > 1 
by using the second line of Eq. (l9|, wc find the local 
diffusion coefficient D{x) ^ e"^'^. It enhances from 

x{L, — x) 

its value in the passive case, which is ^ e ^ [9], 
by a factor of ~ exp[i(y-^=-)^]. The enhancement 
factor increases monotonically from the interface to the 
sample mid-point. Such inhomogeneous enhancement re- 
fiects that upon turning on absorption resonant states are 
'killed', and the portion is determined by the lifetime of 
resonant states. Indeed, near the sample mid-point these 
states have the longest lifetime and are most easily to be 
'killed', and this accounts for the strongest enhancement 
at the mid-point. 

The interplay of localization and absorption is even 
stronger in the latter case {£,a/L <C 1). For x suffi- 
ciently away from the sample center, x <C Ca/2, the 
scaling factor (l8| reduces to that of semi-infinite sam- 
ples, X{x) = 21(1 ~ e-^"^/^"), and D{x) is the same as 
that of semi-infinite samples. (We consider only x < L/2 
since D{x) is symmetric with respect to the sample mid- 
point.) Its behavior is summarized in Table W and ^U 
Near the sample center, x ^ Ca/2, the scaling factor Q 
is simplified to A(a;) w A(L/2). Thus, D{x) w D{L/2) 
exhibiting a plateau. Strikingly, despite that the sys- 
tem is controlled by two parameters, ^/L and ^a/L, the 



TABLE II: Behavior of D{x) at £,a <. ^ <. L 



regime 


D{x)/D(0) 


X > Ca/2 


l+Cl{^a/20+C2{^amr + --- 


X < ^a/2 


l+C^{x/0+C2{x/0'+--- 



plateau value depends on single parameter, ^a/C- More 
precisely, Eqs. (pi) and (Is]) give Eq. pi. Moreover, from 
Eq. ([9| we obtain the expression of D{L/2) for ^a/2^ ^ 1 
(Ca/2C ^ 1) which is given by the first line of Table IT] 
( [lT| . Note that in the strong absorption limit, $,/(,a — >■ oo, 
the homogeneity of the diffusion coefficient is restored, 
D{x) = Dq, reflecting the absence of localization effects. 

That the spatial resolution of D{x) for x < L/2 is 
identical to that of semi-infinite samples when £,a/L <S^ 1 
reflects an important fact: for strong absorption waves 
could not propagate from the left-half to the right-half 
part of the sample and vice versa. As such, resonance 
states play no roles. Therefore, one might expect that the 
phenomenological model of Ref. [TT] is valid, as observed 
in numerical experiments |17) . 

In summary, we present the first microscopic theory 
showing that in localized open absorbing media photons 
display unusual Brownian motion. First of all, the dif- 
fusion coefficient is inhomogeneous in space; most strik- 
ingly, despite that the system is controlled by two pa- 
rameters {£,/L and ^a/L), it exhibits novel single param- 
eter scaling namely Eq. (pj. The analytic predictions 
for the spatial resolution of the diffusion coefficient are 
confirmed by numerical simulations. We stress that our 
theory is very general and, particularly, valid for arbi- 
trary absorption strength. In the limiting case of strong 
absorption realized experimentally [16] . our results agree 
well with experimental measurements. It is very interest- 
ing to generalize the present theory to the gained system 
in the future, and this may have direct applications in 
random lasers [5S]. 
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